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Abstract. Atomic- resolution electron microscope images show that a quasicrystal is a quasiperiodic packing 
of clusters. The outer atomic shells of multi-shell clusters occuring in quasicrystals are highly symmetric and 
rather robust, but some structural disorder may be present around the core. The models describing qua- 
sicrystal structure are quasiperiodic lattices containing two or more unit cells decorated by atomic clusters. 
We show that a modified version of the strip projection method may become an alternate way to obtain 
I . useful models. 
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^ ■ 1 Introduction 

It is known that all quasicrystals found so far are composed of building units, usually called atomic clusters 
Oh" [TJ |2l [31 SI El [HI El [H] - An atomic cluster is a set of close atoms, generally, distributed on fully occupied 

high symmetry shells, which are found at a high frequency in quasicrystal [5J [5] • Atomic-resolution electron 
microscope images offer direct information on the atomic clusters occuring in quasicrystal and cluster packing. 

S The outer atomic shells of multi-shell clusters occuring in quasicrystals are highly symmetric and rather 
robust, but some structural disorder may be present around the core |3]. 
. Any finite set admitting a group G as a symmetry group is a union of orbits of G, and is usually called 

a G-cluster. A mathematical algorithm for generating quasiperiodic point sets by starting from G-clusters 
was proposed by author in collaboration with Verger-Gaugry several years ago |10j . It is based on strip 
projection method and is a direct generalization of the algorithm presented by Katz and Duneau in [llj . 
The model used for the icosahedral quasicrystals in [TT] starts from the one-shell icosahedral cluster C formed 
by the vertices of a regular icosahedron. The physical space is embedded into the superspace R 6 such that 
the orthogonal projections on the physical space of the points 

(±1,0,0,0,0,0), (0, ±1,0, 0,0,0), (0,0, 0,0, ±1,0), (0,0,0,0,0,±1) 

are the vertices of a regular icosahedron. 

In our direct generalization, we consider only G-clusters invariant under inversion. If our starting G- 
cluster C has 2k points then we embed the physical space into the superspace R fe in such a way that C is the 
orthogonal projection on the physical space of the subset 

{(±1,0,0,...,0), (0,±1,0,0,...,0), (0,0,...,0,±1,0), (0,0,...,0,±1)} 

of M fe containing 2k points. Strip projection allows to obtain a pattern V such that for each point x € V 
the neighbours of x belong to the translated copy x + C of C. This means that V is a quasiperiodic packing 
of interpenetrating partially occupied copies of C. In the case of a multi-shell cluster our algorithm uses 
a superspace of rather high dimension and the occupation of clusters occurring in the obtained pattern is 
extremely low. 

Quasicrystals are materials with perfect long-range order, but with no three-dimensional translational 
periodicity. The description of atomic structure of quasicrystals is a highly non-trivial task, and we are 
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still far away from a satisfactory quasicrystal structure solution [J. The quasicrystal structure is usually 
described by using a quasiperiodic lattice containing two or more unit cells decorated with atoms. This 
description is purely mathematical and does not provide any physical insight on why the atoms should 
favour such a complicated structure (3j- Our aim is to present a modified version of the strip projection 
method which leads to useful models and might be justified from a physical point of view. 

2 Packings of G-clusters defined in terms of the strip projection 
method 

Let G = D n be one of the dihedral groups Dg, -Dio, D\ 2 . The group G can be defined as 

G = ( a, b | a n = b 2 = {abf = e ) (1) 

and the formulae 



a(a, 0) = (a cos % - sin ^, a sin % + cos ^) 
b(a,0) = (a,-0) 



(2) 



define an R- irreducible representation in K 2 . Let 

C = {vi, v 2 , v k , -vx, -v 2 , -v k } (3) 

where v\ = (un,W2i), v 2 = {vi2, i>22)r"> v k = (vik, «2fe), be a fixed G-cluster symmetric with respect to the 
origin. From our general theory [TU] (a direct verification is also possible) it follows that the vectors 

w 1 ^ (v n ,v 12 , ...,v lk ) and w 2 = (v 2 i, v 22 , — , v 2 k) (4) 

from M. k are orthogonal and have the same norm 

(wi,w 2 ) = vnv 2 i + v 12 v 22 + ... + vi k v 2k = 

(5) 



IKII = V v ii + v i2 + ■■■ + v tk = \J v i\ + v ii + - + v lk = IKII- 

We identify the physical space with the two-dimensional subspace 

E = { a Wl + 0w 2 | a, S R } (6) 
of the superspace M fe and denote by E 1 - the orthogonal complement (figure 1) 

E-t- = { x eR k \ (x,y) =0 for all y € E }. (7) 
The orthogonal projection on E of a vector x € R fe is the vector 



W! \ Wi I w 2 \ w 2 

X, — > h < X, — 

K I K \ K 



(8) 




Figure 1: Left: The strip S = E + [-1/2, 1/2] 3 and the window W = 7r- L ([-l/2, 1/2] 3 ) in the case of a 
one-dimensional physical space E embedded into the three-dimensional superspace R 3 . Right: The one-shell 
Cg-cluster C = Cg(l, 0) and a fragment of the set V defined by using this cluster and strip projection method 
in a four-dimensional superspace. The nearest neighbours of any point p G V belong to p + C, which is a 
copy of C with the centre at point p. The centres of fully occupied clusters are indicated by o. 
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Figure 2: Left: A fragment of the set defined by using the one-shell Cio-cluster C = Cio(l,0) and strip 
projection method in a five-dimensional superspace. This pattern does not contain fully occupied clusters, 
and its points represent the vertices of a Penrose tiling. The centres of the clusters with occupation greater 
than 50% are indicated by o. Right: A fragment of the set defined by using the one-shell Ci2-cluster 
C = 612(1, 0) and strip projection method in a six-dimensional superspace. The centres of the clusters with 
occupation greater than 50% are indicated by o. 




Figure 3: Left: A fragment of the set V defined by using the two-shell cluster C = Cio(l, 0)UCio(l-l, 1.3) and 
strip projection method in a ten-dimensional superspace. The nearest and the ncxt-to-nearest neighbours of 
any point p € V belong to p + C, but the occupation of the copies of C occurring in this pattern is extremely 
low. Right: The superspace decompositions R k = E © E 1 - = E © E' © E" = £ © E". The points of 7L k lie in 
a family of subspaces parallel to £ — E © E' . 



where k = \ \wi\ \ = \\w2W, and the orthogonal projector corresponding to E is 



E A 



TT X = X — TT X. 



(9) 



We describe E by using the orthogonal basis {k 2 wi, k 2 u>2j-, and therefore, the expression in coordinates 
of 7r is 

7T-.R k >R 2 TTX = ((X,W 1 ),(X,W 2 })- (10) 

The set defined in terms of the strip projection method [111 1121 1131 114] 

v = -k(s n z fc ) = { ttx \ x e s n z k } 

by using the strip (figure 1) 

S = A + E = {x + y\ xeA, y E E } 
generated by shifting along E the unit hypercubc 



A = <^ {xi,x 2 , Xk) 



~^<Xi<- far alii €{1,2,..., fc} 



(11) 
(12) 

(13) 
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is a packing of partially occupied copies of C. The set of arithmetic neighbours of a point x € Z fe is 

A(x) — { x + ei, x + e 2 , £ + e fc , x - ei, a; - e 2 , x - et } (14) 

where ei = (1,0, ...,0), e 2 = (0, 1, 0, 0), = (0, ...,0, 1) are the vectors of the canonical basis ofR fc , 

and the set of neighbours of a point 7rx S is n(S D .A(:r)). Since 

7re 4 = ((e„uii), (e,, iu 2 )) = («it,«2i) = ^ (15) 

we get 

7r(iS (~1 »4(x)) C {7ra: + Wi, 7rx + w 2 7rx + Vk , 7rx — v\ , 7re — w 2 ttx — Ufc} = 7rx + C (16) 

that is, the neighbours of any point ttx G P belong to the translated copy 7rx+C of C. A larger class of 
quasiperiodic patterns can be obtained by translating the strip S. The set 

V = n((t + S) n Z fc ) = { ttx | i-fe5 and i £ Z 1 ' } (17) 

is a packings of partially occupied copies of C, for any ieR fc . 

In the right hand side of figure 1 we present of a fragment a pattern V defined by starting from the 
C$ -cluster 

C=C 8 (1,0) = (±(1,0),± (^j=,^j=) ,±(0,1), ± r ~ l 1 



In this case, the dimension of the superspace is k = 4, 



dimE- 1 = dim_E = 2, and one can remark that V contains fully occupied clusters. The points of V are the 
vertices of an Ammann-Beenker tiling. 

The set from the left hand side of figure 2, obtained by starting from the Cio-cluster C = Ciq(1,0), is 
formed by the vertices of a Penrose tiling and does not contain fully occupied clusters. In this case the 
dimension of the superspace is k = 5 and dimE- 1 = 3=1 + dimi?. If we start from the Ci 2 -cluster 
C = Ci 2 (l, 0) then we obtain a set which does not contain fully occupied clusters (right hand side of figure 
3). In this case k = 6 and dimE^ =4 = 2 + dimE. The set from the left hand side of figure 5 is generated 
by starting from the two-shell Cio-cluster C — Cio(l,0) U Cio(l.l, 1.3). The nearest and the next-to-nearest 
neighbours of any point u> belong to to + C, but the occupation of the copies of C occurring in this pattern 
is extremely low. The dimension of the superspace used in this case is k = 10 and dimi?^ =8 = 6 + dimE 1 . 
These examples show that the occupation of the clusters occurring in our patterns diminish quickly when 
dimE" 1 - — dimE increases. 

One can prove [14] that the space E 1 - can be decomposed into a direct sum E 1 - = E' © E" such that Z k 
is contained in a discrete family of affine spaces parallel to £ = E(&E' (see right hand side of figure 3). Only 
a finite number of these subspaces meet the window W = tt (A) corresponding to the strip S = A + E. In 
the case when dimE^ = dimE the space E" is the null space {0} and the set n (S (1 7L k ) is dense in the 
window W = ^^(A) corresponding to the strip S. This explains why we get fully occupied clusters in the 
case of C = E> 8 (1, 0) and in the case of the model used by Katz and Duneau [11] and independently by Elser 
[T2] for icosahedral quasicrystals. 



3 Modified strip projection method 

The distances to the 'physical' space E of the points of Z fe form a discrete sequence. In table 1 we present 
the smallest terms of this sequence in the cases presented in figures 1 and 2. One can see why the occupation 
of clusters is decreasing when we pass from Cs(l, 0) to Cio(l, 0) and Ci 2 (l, 0). 

The left hand side of figure 3 suggests us that in the case of the sets V defined by multi-shell G-clusters 
there is no cluster occuring at a high frequency. Therefore, generally, the algorithm presented in previous 
section does not lead to useful models. In addition, this purely mathematical algorithm does not provide 
any physical insight on why the atoms should favour such a structure. In this section our aim is to present 
a modified version which leads to interesting quasiperiodic sets, and might be justified from a physical point 
of view. 
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Table 1: Distances to the space E of the nearest points of 7L k in the case of octagonal, decagonal and 
dodecagonal clusters used in figures 1 and 2. 

We start from a G-cluster and embed the 'physical' space E into the superspace R fc in the same way as 
in the previous section. Then we translate E by a certain vector t € M fc and consider the set 

£ = {xeZ k | \\x-t\\<R} 

of all the points of Z k lying inside the ball of centre t and a fixed radius R. For each point x € C we 
determine the distance from x to the space t + E. 

We define a subset Q C n(Z k ) by starting from the assumptions: 

- the distance between any two points of Q is greater than a fixed minimal distance 5 

- each point q = nx from Q has the tendency to become the center of a translated copy of C, and this 
tendency diminishes when d(x) increases. 

We generate a fragment of Q by considering the points x of C in the increasing order of d(x). If the distance 
from ttx to one of the already obtained points of Q is less than S then we pass to the next point of £. If the 
distance from nx to the already obtained points of Q is greater than S then we add irx to Q and continue 
by adding all the points of the cluster ttx + C satisfying the same restriction (distance from the considered 
point to the already obtained points of Q is greater than 5). Then we pass to the next point of £. After a 
certain number of steps the retriction concerning the minimal distance does not allow us to add any point. 

The set Q from figure 4, obtained by starting from the cluster C = £>io(l,0), contains a significant 
percentage of fully symmetric clusters, and the simulated diffraction pattern consists in sharp Bragg peaks 
disposed according to D w symmetry. The minimal distance imposed in this case is the distance between two 
neighbouring points of C. The remarkable set presented in figure 5 is obtained by starting from the cluster 
C = Di 2 (l, 0) by imposing also 5 to be the distance between neighbouring points of C. The starting cluster is 
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Figure 4: Left: Quasiperiodic packing of one-shell decagonal clusters obtained by using the modified strip 
projection method. The centers of the fully symmetric clusters are indicated by o. Right: The simulated 
diffraction pattern. 
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Figure 5: Left: Quasiperiodic packing of two-shell dodecagonal clusters obtained by using the modified 
projection method and the one-shell £>i2-cluster C = -Di 2 (l, 0) as starting cluster. The centers of the clusters 
with outer shell fully occupied are indicated by o. Right: The simulated diffraction pattern contains only 
extremely sharp Bragg peaks, disposed acording to the D±2 symmetry. 

a one-shell cluster but in Q we found a packing of two-shell Z?i2-clusters with the outer shell fully occupied 
and the inner shell containing only 5-6 points. The simulated diffraction pattern corresponding to Q contains 
only extremly sharp Bragg peaks and admits D12 as a symmetry group. The computer programmes used in 
this paper are available on line |15j . 
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